Two Component Heat Diffusion Observed in CMR Manganites 
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We investigate the low-temperature electron, lattice, and spin dynamics of LaMnOa (LMO) and 
Lao.7Cao.3Mn03 (LCMO) by resonant pump-probe reflectance spectroscopy. Probing the high-spin 
d-d transition as a function of time delay and probe energy, we compare the responses of the Mott 
insulator and the double-exchange metal to the photoexcitation. Attempts have previously been 
made to describe the sub-picosecond dynamics of CMR manganites in terms of a phenomenologi- 
cal three temperature model describing the energy transfer between the electron, lattice and spin 
subsystems followed by a comparatively slow exponential decay back to the ground state. However, 
conflicting results have been reported. Here we first show clear evidence of an additional component 
in the long term relaxation due to film-to-substrate heat diffusion and then develop a modified three 
temperature model that gives a consistent account for this feature. We confirm our interpretation by 
using it to deduce the bandgap in LMO. In addition we also model the non-thermal sub- picosecond 
dynamics, giving a full account of all observed transient features both in the insulating LMO and 
the metallic LCMO. 

PACS numbers: 75.47.Lx, 72.15.Eb, 72.25.Rb, 78.47.D- 



I. INTRODUCTION 

Transition-metal oxides with partially filled d-electron 
bands exhibit a wide range of fascinating phenomena like 
dynamic lattice distortions, magnetic, charge, and or- 
bital ordering. Particularly the perovskite manganese ox- 
ides La(!_ x )A x Mn03, where A is a divalent alkali metal, 
have attracted scientific attention. While the undoped 
LaMnOs (LMO) is an A-type antiferromagnetic Mott 
insulator^ 2 - of orthorhombic symmetry, Lao.7Cao.3Mn03 
(LCMO) exhibits the colossal magnetoresistance (CMR) 
effect along with the transition from a paramagnetic in- 
sulating to a ferromagnetic metallic phase 3 -. It is ac- 
companied by a structural transition from orthorhom- 
bic to rhombohedral symmetry. With respect to techni- 
cal applications of the CMR manganites, for instance in 
magneto-optical switches, the dynamics of their strongly 
correlated electron, lattice, and spin subsystems are of 
particular interest^. 

Pump-probe spectroscopy provides a powerful tech- 
nique to investigate these dynamics: A short fs pump- 
laser pulse yields instantaneous photoexcitation of the 
electron system via electric-dipole allowed transitions. 
It has been established that, after the restoration of a 
Fermi distribution by electron-electron interactions (< 1 
ps), excess heat is transferred from the electron system 
to the lattice on a ps time scale via phonon emission 5 . 
Successively, heat is transferred from the lattice to the 
spin system on a ps to ns time scale mediated by the 
spin-orbit interaction- - - 7 -. The dynamics of the electron, 
phonon, and spin systems are probed by a second laser 
pulse delayed with respect to the pump pulse. 

Previous pump-probe studies of the CMR manganites 
have focused on the spin-lattice (1-s) relaxation while 



the dynamics of the electron-electron (e-e) and electron- 
lattice (e-1) interactions have been largely neglected so 
far. For instance, the temperature dependence of the 
transient transmission of thin LSMO and LCMO films 
was analyzed in terms of two components with time scales 
< 1 ps and 20-200 ps. While the sub-ps transient is natu- 
rally associated with e-e processes, the slower component 
was attributed to the 1-s relaxation 8 . 

Recently however, differential transmission spec- 
troscopy studies on the antiferromagnetic parent com- 
pound LMO explained the neglected intermediate tran- 
sients in terms of the relaxation of the collective Jahn- 
Teller distortion and energy transfer between the elec- 
tron and phonon subsystems. Curiously though, the 1-s 
component seemed absent^. Similarly, the correspond- 
ing relaxation in LCMO was attributed to Jahn- Teller 
polarons~ There are, however, features in the picosec- 
ond dynamics in CMR manganites not covered in the 
published literature. One such feature is a sign change 
after several hundreds of picoseconds 11 '. With the spin- 
lattice relaxation time in the range 10-100 ps, such a 
zero crossing is not possible according to the models pre- 
viously used. Also, the observed dependence of the heat 
diffusion rate on the probe energy 7 is in contrast to the 
models used. 

In this paper, we study the free carrier, electron, lat- 
tice, and spin dynamics of LMO and LCMO by reso- 
nant pump-probe reflectance spectroscopy. A modifica- 
tion of the widely used phenomenological three temper- 
ature model first described by Beaurepaire^ is developed 
in order to explain all the observed features, including 
the aforementioned features previously unaccounted for, 
in both LMO and LCMO. For delays larger than 1 ps 
we observe three clearly distinguishable components in 
the low-temperature transients. First we see a fast com- 
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ponent reflecting the transient electron population after 
photoexcitation (~ 5 ps) followed by an intermediate 
component reflecting the lattice temperature dynamics, 
exponentially rising with the e-1 relaxation time, expo- 
nentially decaying with the 1-s relaxation time. Sub- 
sequently we have a slow component (exponential rise 
time ~ 80 ps) in agreement with previous reports of 
conventional 1-s heating by spin-orbit interaction. We 
conclude that the intermediate component relates to the 
presence of exchange mechanisms in the ground state. 
In CMR compounds, lattice heating induces a relaxation 
of the Jahn- Teller distortion which causes a change in 
the transfer matrix of the exchange mechanisms. Fol- 
lowing these relaxation components we observe finally a 
decay on the ns scale consisting of two components. This 
bi-exponential decay back to the ground state has been 
mentioned before^ but without much discussion. Here 
we show a clear separation of the two processes where 
the faster of these is interpreted as the film-to-substrate 
heat diffusion (~ 500ps) while the slower (~ 10 ns) is 
interpreted as the in-film heat diffusion. 

The sub-picosecond response due to changes in the 
density of states immediately following the photoexcita- 
tion is seen as a sharp peak centered around zero delay, 
and the rise and decay of this peak is modeled by a phe- 
nomenological differential equation that almost exactly 
reproduces the observed response. Solving the differen- 
tial equation, the response is basically a convolution be- 
tween the pump pulse and an exponential decay (~ 70fs) 
representing the thermalization of the excited electrons 
back to a well defined fermi level. There is also signs of a 
photoinduced phase transition in LMO where the ~ 70fs 
component of metallic origin is accompanied by a slower 
(~ 850 fs) component indicating the reformation of the 
band gap. Using the three temperature model modi- 
fied with two component heat diffusion together with the 
modeled non-thermal response, all observed features in 
the transients of CMR manganites can be understood. 



II. EXPERIMENTAL 

The transient reflectivity measurement was performed 
at near normal incidence with a Tksapphire mode-locked 
amplifier with a 1kHz repetition rate and nominal 170 
fs pulse length. The pump beam uses the fundamental 
at 775 nm while the probe energy E pr is tuned over the 
range 0.59-2.25 eV using an optical parametric amplifier 
(TOPAS). The beam diameters at the sample surface was 
adjusted to 135 and 160 fim for the probe and pump re- 
spectively, giving a constant pump fluence < lmJ/cm 2 
and a probe/pump fluence ratio of less than 1/30. Us- 
ing the cross-correlation between pump and probe, the 
effective time resolution is estimated to ~ 190 fs, while 
autocorrelation gives an individual pulse length close to 
the nominal 170 fs. 

The LMO and LCMO films of respective thick- 
nesses 175 and 165 nm were epitaxially prepared on 



SrTiO 3 (100) and NdGaO 3 (110) substrates by pulsed- 
laser deposition using stoichiometric targets. The c 
direction of the LMO film is orthogonal to the surface 
normal. In the A- type antiferromagnet LMO, the Mn 
spins order ferromagnetically in the ab plane while they 
are stacked antiferromagnetically along the c direction. 
With the electric-field vector parallel to the ab plane we 
thus excited and probed ferromagnetic spin correlations 
in this experiment. The Neel and Curie temperatures 
of the LMO and LCMO films are T N - 150 K and 
Tc = 256 K, respectively. The metal- insulator transition 
temperature of the latter is Tmi = 264 Ki2. 



III. RESULTS 

The reflectivity transients of LMO between 1.31 and 
2.25 eV at 80K are shown in Fig. QJa). After the initial 
negative peak at t = the transient increases, changes 
sign between ~5 ps (Epr = 2.25 eV) and ^100 ps (1.48 
eV), reaches a maximum between ^400 ps (2.25 eV) and 
beyond 2ns (1.48 eV), and eventually decays exponen- 
tially on a ns time scale. We model the experimental 
results using a modified three temperature model. The 
standard three temperature model^ describes the energy 
flow into and from the electron, lattice and spin subsys- 
tem respectively, followed by a usually very slow heat 
diffusion component: 

dR/R = A e e-*/ Te! +Ai{l- e -*f T A e~ tlru 

+A S {l - e - t/T '^ e- t/Td (1) 

where -A 6) ; jS describes the magnitude on dR/R arising 
from an increased temperature in the electron, lattice 
and spin system respectively. Similarly, Td,ei,u denotes 
the time constants for heat diffusion (ns-scale), energy 
transfer between the electron and lattice systems (ps- 
scale) and lattice and spin systems respectively (100 ps 
scale). A different interpretation for the ns-scale pro- 
cess has been suggested^ where the slowest relaxation 
back to the ground state was attributed to metastable 
quasiparticle dynamics involving spin-flip processes. As 
the slowest component we observe is present also in the 
paramagnetic room temperature phase (not shown here) 
we keep the orthodox heat-diffusion interpretation. The 
described model should be viable as long as only ther- 
mal processes occur. It has been previously reported 
in similar materials that this is the case after ~ lps, 
after which the initial non-equilibrium electron gas has 
had time to thermalize back to a fermi-distribution. It 
should be noted that the assumption of thermal processes 
implies that the time constants for the energy transfers 
between the subsystem as well as the subsequent heat dif- 
fusion will be the same regardless of probing energy. A 
systematic fitting procedure, in which the time constants 
were held constant between all probe energies, were per- 
formed on all our data. It is possible to fit individual 
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FIG. 1: (Color online) Transient reflectivity in (a) LMO and 
(b) LCMO for different probe energies when pumping with 1.6 
eV. Notice in particular the change in the final relaxation time 
when changing probe energy and the zero crossing (~ 100 
ps, LMO and ~ 550 ps LCMO) when probing with 1.47eV, 
features that are so far unexplained by the literature. 



transients using models with fewer components, but this 
simplification comes with the expense that the thermal 
interpretation with time constants independent of prob- 
ing energy has to be abandoned. With this systematic 
fitting procedure it is possible to fit all LMO data using 
an e-1 time constant r e ; = 10 ps and a 1-s time constant 
Ti s = 75 ps. Our value of r e ; is in good agreement with 
previous differential transmission studies, though no sign 
of the spin dynamics that we link to the A s component 
was reported^. In addition, a consistent fit to our data re- 
quires two slow components t<h = 720 ps and T42 = 8960 
ps: 



dR/R = A e e- t/T " 1 + Ai ( 1 - e"* /T ' 



1 ) p-t/TU 



+A S (l - e-'/^) e- t/Tdl + C r e~ t/Tt!2 (2) 

where C r is the relative weight between the two slow 
components. The necessity of both these components 




FIG. 2: (Color online) The computed fittings to the experi- 
mental (a) LMO and (b) LCMO data plotted on a semi-log 
scale. The dots represent the experimental data while the 
lines is the corresponding fittings. The inset show the same 
fittings using a linear scale. The transients shown are the 
same as in Fig. [T] 



can be directly deduced from Fig. QJa) by looking at the 
1.3 and 1.47 eV data. First of all, the differential reflec- 
tivity seen by the 1.3 eV probe shows a, in comparison, 
very fast decay to zero signal, clearly not compatible with 
what is seen at higher probe energies. Also, the 1.47 eV 
transient has its maximum after a delay larger than 1500 
ps, indicating at least two components at interplay after 
several hundreds of picoseconds. 

The corresponding reflectivity transients in Ca-doped 
LCMO at 80 K are displayed in Fig. [jib). Apart from 
E pr = 1.65 and 1.47 eV, where the sign changes around 
0.6 (see Fig. HJb)) and 550 ps, respectively , the ini- 
tial peak and long-lived transient show the same sign 
depending on E pr . Note that above t ~ 1 ps, the tran- 
sient increases or decreases much faster than in LMO, 
reaching an extremum between about 70 ps (2.25 and 
1.31 eV) and roughly 400 ps (1.65eV) before exponen- 
tially decaying within several ns. Additionally, the tran- 
sient is modulated by a damped oscillation below ~ 40 
ps. The same modified three-temperature model intro- 
duced above, Eq. is sufficient also to describe the 
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low-temperature transient reflectivity of LCMO, with the 
addition of an oscillating part attributed to the release of 
coherent acoustic phonons^. The resulting fit is shown 
in Fig. [5Jb) and gives for LCMO the time constants 
r el = 2.6 ps, t 1s = 31.5 ps, r dl = 475 ps, r d2 = 7000 
ps and an oscillation period E pr T p ~ 60eV ps. This in- 
dicates a faster response than in LMO, which may be 
attributable to the stronger correlations in LCMO. Our 
observed time constants are consistent with the values 
obtained from time resolved conductivity measurement!^ 
and time resolved MOKE measurements^ 7 -^, although 
strain effects due to the use of different substrates might 
influence the exact numbers. 

In Fig. [3]we show the values of the magnitudes A e j tS ^ 
and C r obtained for LMO and LCMO from the fitting 
procedure described above. A e , Ai and A s all show very 
similar trends with probe energy, while the dependence 
of C r on probe energy is very different. The accuracy of 
the parameters is best for A s and C r due to the large 
signals and decoupling in time from the other processes. 
In LCMO, accurate determination of A e and Ai is hin- 
dered by the simultaneous presence of oscillations from 
the coherent phonons, something that is reflected in the 
large errorbars for these parameters. In both LMO and 
LCMO, C r shows a divergent behaviour. The divergence 
occurs when approaching the energies 1.2eV (LMO) and 
1.4eV (LCMO) and is discussed below in Section llVl for 
the insulating LMO, while so far the divergence in LCMO 
is unexplained. 

Having successfully modeled the transients above the 
picosecond scale with the three temperature model, we 
turn now to the non-thermal sub-picosecond response. 
Typical ultrafast transients for LMO and LCMO is shown 
in Fig. [4] (a) and (b) respectively. The main shape of the 
peak is gaussian with a width corresponding to the cross- 
correlation between the pump and probe pulses. Con- 
voluted with the gaussian shape is a fast decay due to 
the thermalization of the excited electrons, while LCMO 
also shows high frequency oscillations with a period 
~ 600fs, possibly from strongly damped coherent opti- 
cal phonons 17 . To account for these features, neglecting 
the oscillations, we model the number of electrons N in 
an excited state as 
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FIG. 3: (Color online) The computed fitting parameters of 
LMO (blue circles, right axis) and LCMO (red triangles, left 
axis) where the errorbars indicates 95% confidence intervals. 



dN 
dt 



p(t) - N/r t 



(3) 



where t± is the decay rate (thermalization time) and 
P(t) ~ cxp (— i 2 /Tp) is the source term with a width t p . 
Solving Eq. © and assuming further that dR/R ~ N, 
the non-thermal contribution to the transient response 
gets the form 

dR/R = i [1 - erf (r p /2r t - t/r p )} A t e^ Tt (4) 

with erf the error function and A t a parameter to be 
determined, together with r t , by fitting to the data. Re- 
quiring the fittings in Eq. @ and (@| to be consistent, 



a single equation describing the complete transient re- 
sponse is obtained 



dR/R = i [1 - erf (r p /2r t - t/r p 
A t e' t/Tt + A e (l - e~*/ T *) e- t/T " + 



(5) 



where . . . denotes the remaining part of Eq. ([2]). 

When fitting the peak and subsequent carrier thermal- 
ization the duration t p of the pump pulse is assumed to 
be the nominal 170 fs. Using the parameters obtained 
from Eq. @, there are two parameters, Tt and A t , to 
be determined by the same fitting procedure as above. 
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FIG. 4: (Color online) Sub-picosecond transients in LMO (a) 
and LCMO (b). In both samples we see a roughly gaussian 
peak followed by very fast (< 1 ps) decay processes, while 
only LCMO shows signs of coherent optical phonons. The 
blue dots represents the experimental data while the red lines 
are the corresponding fittings. 



Using this scheme, the thermalization time Tj is found 
to be 67fs for LMO and 78fs for LCMO. The amplitude 
At, describing the impact of the excitations of the free 
carriers, exhibits an energy dependence similar to that of 
A e , Ai and A s . 

The response modeled in Eq. ([5]) is completely ade- 
quate for LCMO but for LMO an additional fast compo- 
nent with decay time 850 fs is needed for a satisfactory 
fitting. The existence of a two component decay of pho- 
toexcited carriers in Mott insulators have been reported 
previously^ and has been attributed to an immediate re- 
sponse following the photo induced insulator-metal tran- 
sition followed by a transition to a slower relaxation after 
the recovery of the insulating bandgap. 



IV. DISCUSSION 

In pump-probe spectroscopy, electric-dipole allowed 
transitions are excited and probed. As their strength 



depends on the densities of initial and final states 
(DOS) and the transition probabilities between them, 
reflectivity transients arise from either DOS changes or 
transition-probability changes. From numerous stud- 
ies of the equilibrium optical propertie a 2 ' 19 ' 20 and band- 
structure calculations 21 of the manganites two kinds of 
electric-dipole allowed transitions may be expected in the 
spectral range 0.8 — 2.3 eV: transitions between the 3d 
orbitals of neighboring Mn ions (d-d) and transitions be- 
tween 02p and Mn3d states (p-d). While the former are 
highly sensitive to the relative orientation of Mn spins 
and thus on temperature, the latter are expected to be 
temperature independent^! . 

While the initial photoexcitation by the pump pulse 
might involve both d-d and p-d transitions, the strong 
temperature dependence of the intermediate and slow 
transient components, especially across Tn and Tc, 
clearly indicates that these components arise from prob- 
ing d-d transitions£i£i&2£. 

Let us summarize the current understanding of pho- 
toinduced dynamics in LMO or LCMO. The 1.6 eV pho- 
toexcitation creates pairs of electrons and holes (quasi- 
particles). Within < lps, the electrons excited into the 
upper d band thermalise to a Fermi distribution at an ele- 
vated temperature T e , out of equilibrium with the lattice 
and the spin system, occupying formerly empty states at 
the lower band edge ». We found that this process is a two 
component process in the insulating LMO (r = 67 and 
850 fs respectively), while a single component is sufficient 
for LCMO (t = 78)fs. Subsequently, the quasiparticles 
recombine via broad-band phonon emission (e-1 relax- 
ation). The lattice temperature rise at the expense of the 
electron temperature. In the next step, any excess lattice 
energy is transferred to the spin system (1-s relaxation) 
via spin-orbit interaction, the lattice and spin temper- 
atures equilibrate^. This effect was observed in CMR 
manganites before, the s-1 relaxation time was found to 
be ~ 20 to 200 ps, depending on temperature^. 

At this point, the temperatures of all three baths are 
in equilibrium. The temperature in the illuminated vol- 
ume decreases via heat diffusion. We showed above that 
this heat diffusion consist of (at least) two components. 
To interpret these two processes we turn to LMO. The 
relevant band in LMO that is probed and pumped in 
the visible range is the high spin e g — e g charge transfer 
between adjacent Mn 3+ ions via the bridging O 2- ion. 
This charge transfer is mediated both by super-exchange 
and double-exchange, hence being very sensitive to the 
Mn— O— Mn bond angle, indicating the importance of the 
static Jahn- Teller distortion. Thus the spin-temperature 
dynamics become visible because the probed high spin d- 
d transition acts as a thermometer: spin disorder results 
in a weakening of the transition probability, while the lat- 
tice dynamics become visible because an increased lattice 
temperature causes a relaxation of the Jahn- Teller distor- 
tion which thereby changes the Mn — O — Mn bond angle. 
The optical absorption in this band decreases abruptly 



around 1.5 eV— , i.e with E t 



1.3 eV LMO is more 
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FIG. 5: (Color online) The measured coefficient C T compared 
to the theoretical curve (E — Eo)~ 2 on a log-log scale. 



transparent, indicating that the film-substrate interface 
is more important. 

Hence we attribute the fast heat diffusion process to 
the film-to-substrate diffusion, a conclusion shared with 
reference^. Analogously, with higher E pr LMO is more 
opaque, giving us reason to attribute the very slow re- 
laxation process to the traditional in-film heat diffusion 
in accordance with previous assumptions^^. 

To further motivate this interpretation, consider the 
coefficient C r . The magnitude of this coefficient should 
relate inversely to the number of high spin d-d excita- 
tions probed at a given energy. As LMO is an indirect 
semiconductor—, the probability of such an excitation, 
assuming parabolic density of states close to the opti- 
cal gap, should be proportional to (E — E ) 2 where E 
is the probe energy and E$ is the energy of the band 
gap plus the mediating phonon^. Therefore we expect 



C r ~ (E — E )~ 2 . In Fig. [5] we have plotted C r as a 
function of energy together with (E — Eq)~ 2 - Eq was 
obtained by fitting the data to the model, resulting in 
Eq = 1.195eV. This value should be compared to the 
previously published value E gap ~ 1.2eV— . The good 
agreement between the value obtained here with the pre- 
viously published value, together with the good quality 
of the fit, shows that the film-to-substrate interpretation 
is reasonable. 

A similar analysis can be made for LCMO, but be- 
cause LCMO is metallic at 80K, the lack of a band-gap 
never causes us to expect a diminishmcnt in relative in- 
tensity of the in-film heat diffusion when lowering E pr . 
Hence, the film-to-substrate component will only be ex- 
plicit near a sign change, as is the case for E pr = 1.47 
eV, see Fig. [IJb). Keeping the above discussion in mind, 
we can understand the absence of LMO spin dynamics 
in previous differential transmission studies^ as due to 
the thickness of the thin films used there (70 nm), thus 
increasing the effect of the film-to-substrate diffusion for 
all probe wavelengths. 

In summary, we have performed resonance pump- 
probe studies on the femtosecond scale on both the un- 
doped LaMn03 and the doped Lao.7Cao.3Mn03. In spite 
of having very different properties, with LMO being an 
insulator and LCMO being a metal at 80K, the observed 
differential reflectivity transients are very similar in char- 
acter. We have developed one general mathematical 
model, Eq. (O, which captures all observed phenomena 
in the transients, including free carrier, electron, lattice 
and spin dynamics. This model differs from the previ- 
ously used in that it includes two different heat diffusion 
channels, in-film diffusion and film-to-substrate diffusion. 
We conclude that the penetration depth is generally large 
enough for the film-substrate interface to have an effect, 
and that this effect can be used to deduce the bandgap 
in semiconducting samples. 
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